Electron-Phonon Interaction and Ultrasonic Attenuation in the Ruthenate and 

Cuprate superconductors 
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This article derives an electron-phonon interaction suitable for interpreting ultrasonic attenuation 
measurements in the ruthenate and cuprate superconductors. The huge anisotropy found experi- 
mentally (Lupien et at, 2001) in S^RuCU in the normal state is accounted for in terms of the layered 
square-lattice structure of S^RuCU, and the dominant contribution to the attenuation in S^RuCU 
is found to be due to electrons in the 7 band. The experimental data in the superconducting state 
is found to be inconsistent with vertical lines nodes in the gap in either (100) or (110) planes. Also, 
a general method, based on the use of symmetry, is developed to allow for the analysis of ultrasonic 
attenuation experiments in superconductors in which the electronic band structure is complicated 
or not known. Our results, both for the normal-state anisotropy, and relating to the positions of 
the gap nodes in the superconducting state, are different from those obtained from analyses using a 
more traditional model for the electron-phonon interaction in terms of an isotropic electron stress 
tensor. Also, a brief discussion of the ultrasonic attenuation in UPt3 is given. 

PACS numbers: 74.70.Pd, 74.20.Rp, 74.25.Ld 



I. INTRODUCTION 

This article describes a theory of ultrasonic attenua- 
tion in the ruthenate and cuprate superconductors. It 
has been stimulated largely by the work of Lupien et 
allil which presented the results of detailed experimental 
measurements of ultrasonic attenuation in both the nor- 
mal and superconducting states of Sr 2 Ru04. The goal of 
that work was to use ultrasonic attenuation as a tool to 
gain information on the presence and location of nodes in 
the superconducting gap. Quite unexpectedly, however, 
they found a huge anisotropy in the measured attenua- 
tion, even in the normal state. Some of their results are 
shown in Fig. [j]. Notice that the normal-state attenua- 
tion of the longitudinal wave propagating along the [110] 
direction is lower by a factor of approximately 30 than 
that of the [100] longitudinal wave. Furthermore, the 
attenuation of the transverse wave along the [100] direc- 
tion is lower than that of the transverse wave along the 
[110] direction by a factor of more than 1000. Lupien et 
al. note that their results indicate the need for a new 
theory of the electron-phonon interaction allowing for a 
significant variation for the different sound wave modes. 
They also suggest that the lack of a reliable theory of 
ultrasonic attenuation in Sr2Ru04 has greatly hindered 
the use of this technique as a tool for gaining information 
about the location of the gap nodes in this material. 

The use of ultrasonic attenuation as a tool to locate 
the positions of nodes in the energy gaps of supercon- 
ductors has been described, for example, by Moreno and 
Colemaro. In their model, the phonon strain field cou- 
ples to a stress tensor describing the flow of electron mo- 
mentum. The electron stress tensor is taken to be that 
appropriate to an isotropic electron fluid. This approach, 
which can be traced back to early work on ultrasonic 
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FIG. 1. Experimental data on the mode viscosity for the 
four in-plane sound wave modes taken from Ref. The mode 
viscosity 77 plotted here is related to the attenuation a by the 
formula r\ — apv\ /(2nv) 2 , where p is the density, v s is the 
sound velocity and v is the frequency. All sound wave modes 
have both the direction of propagation and the polarization 
lying in the basal plane. The nature of the mode (L for lon- 
gitudinal and T for transverse) as well as the Miller indices 
of the propagation direction are shown in the figure. 
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attenuation in heavy-ikrmion superconductorscTEl, in s- 
wave superconductorsLrp and in normal metalsEI, has also 
been used in the most recent attempts to understand 
ultrasonic attenuation in S^RuOjilM. The use of an 
electron stress tensor appropriate to an isotropic fluid 
gives an elegant and simple formulation of the theory. 
However, it seems clear that a new approach is needed 
if one wishes to be able to account in detail for the ul- 
trasonic attenuation observed in the recent experiments 
on Sr 2 Ru04. The approach developed in this article 
takes account of the crystalline and electronic structure 
of Sr2Ru04, but also can be shown to reduce to the tra- 
ditional isotropic electron-stress-tensor approach in an 
appropriate limit (this limit not being appropriate for 
Sr 2 Ru0 4 ). 

This article shows that the strong anisotropy of the ul- 
trasonic attenuation in Sr2Ru04 is intimately connected 
with its layered square-lattice structure. We start from 
the idea that the electronic structure of Sr 2 Ru04 in the 
neighborhood of the Fermi surface can be well described 
in terms of a simple tight-binding Hamiltonian,oEj in 
which a central role is played by hopping matrix elements 
describing the hopping of an electron from one ruthenium 
ion to a nearest-neighbor or next-nearest-neighbor ruthe- 
nium ion. To develop a model for the electron-phonon 
interaction, we assume that the hopping matrix element 
depends on the distance between the ions, so that if a 
sound wave stretches the distance between two ions, the 
matrix element describing the hopping of an electron be- 
tween these two ions changes. It is easily seen that a 
transversely polarized sound wave travelling in the [100] 
direction in the ruthenium-ion square lattice does not 
stretch the nearest-neighbor bond between two ruthe- 
nium ions, and hence is not coupled to the electrons (at 
least by a nearest-neighbor coupling). This is the reason 
for the extremely low attenuation of the T100 sound wave 
reported in Ref . [l] (see Fig. |l|) . It is also easily seen that 
for a two- or three-dimensional hexagonal lattice, or for 
a three-dimensional body-centered or face-centered lat- 
tice, or if the next-nearest-neighbor interaction is large, 
this argument does not apply. Thus the highly unusual 
strongly anisotropic ultrasonic attenuation observed in 
Sr 2 Ru0 4 is directly related to the layered square-lattice 
structure of this material. The anisotropy in the attenu- 
ation of longitudinal waves (e.g. for the L100 and LI 10 
waves of Fig. [I]) also finds an elementary qualitative and 
quantitative explanation in terms of the model developed 
in detail below: this explanation depends on the details 
of the Fermi surface geometry and will be given later. 

After developing a detailed model for the electron- 
phonon interaction in Sr 2 Ru04, and confirming that the 
model accounts well for the ultrasonic attenuation in the 
normal state, we proceed to an analysis of the ultrasonic 
attenuation in the superconducting state with the objec- 
tive of gaining information on the positions of nodes in 
the superconducting energy gap. The basic ideas here 
have been clearly set out in the article by Moreno and 
Colemana. They show that for a given sound wave prop- 



agation direction and polarization, the nodes in the gaps 
can be described as either "active" or "inactive" , and 
that the temperature dependence of the ultrasonic at- 
tenuation is very different in the two cases. These ideas 
are described in greater detail below, and exploited to 
gain information on the positions of the gap nodes in 
Sr 2 Ru0 4 . Our calculations are done in the hydrody- 
namic limit. The opposite ("quantum") limit relevant 
for either very pure samples, or very high frequencies, 
has been considered in Ref. [14| 

A consequence of our work is the recognition that 
the traditional isotropic electron-stress-tensor model of 
the electron-phonon interaction makes certain nodes "ac- 
cidentally" inactive for the case of longitudinal sound 
waves. This accidental inactivity is a consequence of the 
fact that the isotropic electron-stress-tensor model of the 
electron-phonon interaction is not sufficiently general to 
describe realistically the electron-phonon interaction. In 
a more realistic model, all nodes are active for longitu- 
dinal phonons. Thus, the isotropic electron-stress-tensor 
model for the electron-phonon model gives misleading 
results for the temperature dependence of the ultrasonic 
attenuation associated with longitudinal phonon modes 
in some cases, and should not be used in studies aimed 
at determining the positions of the nodes in the gap of 
an unconventional superconductor. 

Although the emphasis in this article will be on ultra- 
sonic attenuation in Sr 2 Ru04 because of, the availabil- 
ity of experimental data for this material!], it should be 
noted that our detailed results for the 7 band are also ap- 
plicable (except for the interlayer interaction) to cuprate 
superconductors such as YBa 2 Cu30g+a;; furthermore it 
is a simple matter to develop an appropriate interlayer 
electron-phonon interaction for the YBa 2 Cu30g+a; struc- 
ture analogous to the result for Sr 2 RuC-4 developed here. 
The analysis in Sec. |y| allows us to make some gen- 
eral statements about the symmetry-imposed properties 
of the electron-phonon interaction and their manifesta- 
tions for the ultrasonic attenuation in unconventional su- 
perconductors, in-particular, in UPt3. 

The discoveryEj of superconductivity— in the layered 
perovskite Sr 2 Ru0 4 , and the proposalij that the su- 
perconducting Cooper pairs in that material formed in 
a spin-triplet state has stimulated considerable inter- 
est and study (see Ref. [l?] for a review of the prop- 
erties of Sr 2 Ru04, and Ref. |TJ for a review of the 
symmetry classification and the physical properties of 
unconventional superconductors in general). Until re- 
cently, it was thought that the superconducting order 
in Sr 2 Ru04 could b e . de scribed by the order parameter 
d(k) oc z(k x + ik v )iM3. Because the Fermi surface pL 
Sr 2 Ru04 has a quasi-two-dimensional cylindrical formE£l 
with no points at k x — k y = 0, the Fermi surface for this 
order parameter is fully gapped. More recently, however, 
power-law temperature dependence more characteristic 
of a gap having line nodes have been found in a number of 
experiments, including specific heatcpEiL NQRc3, pene- 
tration depthed, thermal conduct ivitye-al23 and ultrasonic 
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attenuatiorfl. It is not easy to reconcile the presence of 
triplet Cooper pairs, a broken-time-reversal symmetry, 
and line nodes in the gap. A number of proposals to 
do so have nevertheless been made. These include spin- 
triplet states characterized by vectors d(k) of the form 
d(k) oc z(sm(k x a) + ism(k y a)) c3cB and having verti- 
cal line nodes in (100) planes, so-called /-wave states 
characterized by d(klet| r 2| r fcj; + ik y )k x k y or by d(k) oc 
z(k x + ik y )(k x — ky) having vertical line nodes 

in (100) and (110) planes, respectively, by^i-wave states 
characterized by d(k) oc zk z (k x + ik y ) 2 ElfE3 and hav- 
ing horizontal lines nodes in the plane k z = £L by states 
characterized by d(k) oc i(k x + iky) cos(fc z c) cB and hav- 
ing horizontal line nodes in the plane k z = 7r/(2c), or by 
states characterized by d(k) oc z(siuQhra/2) cos(k y a/2) + 
icos(k x a/2)sin(k y a/2))cos(k z c/2) Ej£3 and having hor- 
izontal line nodes in the plane k z — -k/c. The analy- 
sis of the ultrasonic attenuation experiments given below 
allows many of these possibilities to be ruled out, and 
suggests that attention be focused on those possibilities 
characterized by the existence of horizontal line nodes. 

The structure of the article is as follows. Section |l| de- 
velops a tight-binding model of the electron-phonon in- 
teraction accounting for the details of the layered square- 
lattice structure occuring in the ruthenate and cuprate 
superconductors. Section [II evaluates a formula giving 
the ultrasonic attenuation in the normal state of Sr2Ru04 
in terms of the model electron-phonon interaction given 
in Section [fi|, and shows that the extremely strong and 
unusual anisotropy of the attenuation is accounted for by 
the model. Section |y| makes use of the model electron- 
phonon interaction to determine the activity or inactivity 
of the gap nodes for various proposed superconducting 
gap structures for Sr2Ru04 thus allowing statements to 
be made about which of the various proposals for the 
gap structure are consistent with the experimental ultra- 
sonic attenuation data. Section shows that a detailed 
model of the electron-phonon interaction is not neces- 
sary to determine which nodes in the superconducting 
state are active or inactive, by showing how to obtain 
such information from symmetry arguments only. Such 
arguments are sufficiently powerful to be applicable to 
cases where the development of a detailed model of the 
electron-phonon interaction is not available. Our detailed 
model for Sr 2 Ru04 is shown to be consistent with these 
arguments (although the isotropic electron stress tensor 
model is not), and furthermore, some new results re- 
lating to the interpretation of ultrasonic attenuation in 
UPt3 are presented. The Appendix gives a discussion 
of the universality of the low-temperature, temperature- 
independent contribution to the attenuation in the super- 
conducting state, showing that the presence (or absence) 
of universal behavior is associated with the activity (or 
inactivity) of the gap nodes. 



II. ELECTRON-PHONON INTERACTION IN 
LAYERED CUPRATE AND RUTHENATE 
SUPERCONDUCTORS 

The approximate two-dimensional nature@ of the 
Fermi surface of S^RuCU suggests that it can be de- 
scribed, to a first approximation, in terms of electrons 
interacting principally through intraplanar interactions. 
The Fermi surface consists of three sheets (see Fig. ||), 
which can be thought of as being deriYiPrj|fi;om the three 



ruthenium orbitals, d xy , d x 



and dyj 1 



I The Hamil- 



tonian describing, the band structure of a single plane can 
thus be writtentj 
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where c) v n a creates an electron with spin a in the z/th 
(y = {xy, xz,yz}) ruthenium orbital on the ruthenium 
ion at the site n of the ruthenium square lattice. 




FIG. 2. Schematic of the Fermi surface of Sr2Ru2 04 show- 
ing the a, f3 and 7 bands. Also shown are portions of the a 
and P bands that are predominantly xz in character. 

Because of the a z plane of reflection symmetry at the 
center of the RUO4 plane, there is no overlap between 
the xy orbitals and those of xz and yz symmetry. This 
means that one sheet of the Fermi surface (the 7 sheet) 
can be attributed to the xy orbitals, while the a and f3 
sheets result from the hybridization of the xz and yz or- 
bitals. Also, from symmetry, there is no nearest-neighbor 
overlap between the xz and yz orbitals, while the largest 
overlap integral for an xz orbital is expected to be with 
nearest-neighbor ions in the ±x directions, since the lobes 
of these orbitals point in these directions rather than in 
the ±y directions. This means that the xz and yz or- 
bitals form approximately independent one-dimensional 
bands, with the hybridization of these band giving rela- 
tively small perturbations to the energies, except where 
the bands cross. A schematic view of the Fermi surface 
is shown in Fig. |^. 

To derive an expression for the electron-phonon inter- 
action, we assume that the hopping matrix elements in 
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Eq. ([[]) in the vibrating lattice depend on the instanta- 
neous positions of the ruthenium ions, r n — r|i ^ -|- u n , 
where is the equilibrium position of the ion at lat- 
tice site n, and u n is its displacement from equilibrium. 
The lowest order contribution to the electron-phonon 
interaction is found by expanding the hopping matrix 
elements in powers of the ionic displacements u n and 
keeping only linear terms. First consider doing this for 
only the nearest-neighbor hopping matrix elements for 
the xz orbitals that lie in the ±x directions relative to 
each other, since this is a relatively simple effectively one- 
dimensional case. This results in the interaction 

gXZ 

X ( c xz,n,o- c xz,n+a.,o- + c xz ,n+a.,tT C xz,n,a) ■ (2) 

Here is the a:th component of displacement of the 
nth ion in the plane, and u n + a>x refers to the ion one 
Bravais-lattice vector a in the positive x direction from 
the nth ion. Notice that, for reasons of symmetry, this 
expression does not contain displacements normal to the 
bond axis. Thus (even though we have not assumed that 
the forces between the ions are central), only displace- 
ments that stretch the bond distance contribute to the 
electron-phonon interaction in this case. Introducing the 
lattice Fourier transforms of the site variables in Eq. (j^) , 
and summing this Hamiltonian over all planes in the crys- 
tal, gives 



H. 



e—ph 
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±q,qj"^qj C l<,k+q+G,er C ^k,o-, 



q,J,k,o-,G 



(3) 



xz, G is a reciprocal lattice vector, A q j — 
".qj destroys a phonon corresponding to 



where v 
a -q,j + a qji 

wave vector q and polarization j, and N is the total 
number of ruthenium ions in the crystal. The wave vec- 
tors q and k are three dimensional wave vectors and this 
Hamiltonian is for the whole crystal. Here 



q,qj 



<? qj ^(k,q), 



with 



fq,: 



-V2i 



Mv) 



F/(k, q ) = g v Yj& ■ R)Ca • *M)) cos ( k • R )> 



(4) 

(5) 
(6) 



R 



where w qj is a phonon frequency, Vj is the sound velocity 
for a phonon of polarization j, ej(q) is a unit vector in 
the direction of the phonon polarization, M is the total 
mass of the ions in a primitive unit cell, a hat (as in q) 
indicates a unit vector. Since we are interested in this 
article only in low frequency phonons, the expression (jq) 



is given to the lowest non trivial order in the phonon wave 
vector q. The result of this section for the interaction of 
phonons with the xz electrons as described by Eq. (||) , is 
given by Eqs. (||), (||), (||), and (||) with v = xz and with 
the sum over R containing a single term with R equal to 
the Bravais lattice vector a. 

Similarly, the interaction of phonons with the yz elec- 
trons is described by Eqs. (|), (f|), (§), and @ with 
v = yz and with the sum over R containing a single 
term with R equal to the Bravais lattice vector b. Also, 
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yz 



the coupling constant is called g 

The nearest-neighbor interaction of phonons with the 
xy (7 band) electrons is described by Eqs. (||), (Q), (g), 
and (||) with v = xy and with the sum over R containing 
two terms with R equal to the Bravais lattice vectors a 
and b. The coupling constant is called cp ' . 

The next-nearest-neighbor interaction of phonons with 
the xy (7 band) electrons is described by Eqs. (||), (Q), 
(H), and ([|) with v = xy and with the sum over R con- 
taining two terms with R equal to the Bravais lattice 
vectors a + b and a — b. The coupling constant is called 
.</-• 

So far, we have considered only electron-phonon inter- 
actions associated with the stretching of bonds between 
ions lying in a single plane. It is reasonable to expect 
these interactions to be larger than the interplanar ones 
because the separation between RUO4 planes is relatively 
large, although this should be confirmed by comparison 
with experiment. However, the intraplanar interactions 
that we have considered so far do not affect transverse 
phonon modes that have their wave vector in the basal 
plane and their polarization perpendicular to this plane. 
This is because, to first order in the displacements, such 
phonon modes do not stretch the bond distances, and 
more rigorously, because the basal plane is a plane of 
reflection symmetry. Hence, to account for the attenua- 
tion of these phonon modes, it is necessary to consider 
interplanar interactions. 

In describing the interplanar interactions it is reason- 
able to consider first interactions involving xz and yz 
orbitals since these orbitals have lobes sticking out of 
the plane, whereas the lobes of the xy orbitals lie in 
the plane. The fact that the parts of the Fermi sur- 
face that are made up fmrn the xz and yz orbitals show 
the largest corrugations^! along the c-axis direction sup- 
port this consideration.EEl In addition we choose to con- 
sider interactions between ruthenium ions at the cor- 
ner and the body center of the unit cell (which are the 
nearest- neighbor interplane pairs). Unfortunately, the 
xz and yz orbitals are not a convenient basis for de- 
scribing this interaction. For this reason, the new basis, 



c^z = (c xz + c yz )/y/2 and c V z = (- 



t )/y/2 will 



be introduced. Clearly, the £z orbitals will interact well 
with themselves along body diagonals in the directions 
of the vectors a + b + c and —a b + c, while the r\z 
orbitals will interact well with each other along the body 
diagonals —a + b + c and a b + c. 

The nearest-neighbor body-diagonal interactions be- 
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tween £z orbitals are described by Eqs. fl|), (|J), (||), 
and (||) with f = £z and with the sum over R contain- 
ing two terms with d equal to the Bravais lattice vec- 
tors Tj(a + b + c) and i(— a — b + c). Also, the nearest- 
neighbor body-diagonal interactions between rjz orbitals 
are described by Eqs. (||), (Q), (||), and (^|) with v = i]z 
and with the sum over R containing two terms with R 
equal to the Bravais lattice vectors |(— a + b + c) and 
i (a — b + c) . The coupling constants for these two in- 
teractions satisfy g^ z — g r,z . Because the best simple 
approximation to the band structure is in terms of the 
xz and yz states, and not in terms of the £z and rjz states, 
the electron-phonon interaction derived in the £z and r\z 
representation should be transformed back to the xz and 
yz representation. 



III. ULTRASONIC ATTENUATION IN THE 
NORMAL STATE OF SR2RUO4 

As noted in the Introduction, there is an extremely 
strong anisotropy of the ultrasonic attenuation in the 
normal state of S^RuO^. This section shows how this 
anisotropy can be accounted for in terms of the electron- 
phonon interaction just described, and the known Fermi- 
surface geometry of Sr 2 Ru04. 

The attenuation constant aj (q) for an acoustic phonon 
of wave vector q and polarization j is given by ctj (q) = 
{vjTqj) , where Tqj is the phonon lifetime. Because 
the formula for the ultrasonic attenuation in the super- 
conducting state is closely related to that for the normal 
state, the latter more general formula for the supercon- 
ducting state will be given here, and the formula for the 
attenuation in the normal state immediately follows. The 
result of evaluating the phonon lifetime in the hydrody- 
namic limit (the electron quasiparticle mean free path 
much shorter than the phonon wavelength) from the self- 
energy of phonon Green's function in the superconduct- 
ing state is 



1 



fro? 



de 



de 



< ( Ff (k, q)r k Re ^e 2 - |A k | 2 ) . 

J I FS 



(7) 



where Nf is the density of states at the Fermi level in 
the normal state, /(e) is the Fermi distribution function, 
and the Fermi surface average is defined by 



<^ 2 (k,q)> 



FS 



J F^k iq )dS u /v k 
J dSk/vk 



(8) 



The integration over dS^ in this equation is over all (the 
7, the xz and the yz) sheets of the Fermi surface, with 
the electron-phonon matrix element chosen appropriately 
for each sheet. The expression (Q) is valid for all sin- 
glet superconducting phases, as well as for unitary triplet 
phases, for which |A k | 2 = |d(k)| 2 . This formula is similar 



to that employed in Ref. |[ except that in our expression 
the electron-phonon matrix element Fj(k, q) replaces the 
isotropric electron stress tensor of Ref. |[ Also, the for- 
mula has been generalized to be applicable to anisotropic 
multi-sheet Fermi surfaces. The quantity r k is the Bo- 
goliubov quasiparticle lifetime, and for Eq. (Q) to be 
valid, the condition ksT Ti/t^ must be satisfied (see 
Appendix for a more detailed discussion). 

One aspect of the above discussion that is unsatisfac- 
tory is the failure to give a full treatment of the Coulomb 
interaction so that charge neutrality is preserved in the 
distorted lattice that occurs in the presence of a lon- 
gitudinal sound wave. A microscopic treatment of this 
question for our tight-binding multi-band model is be- 
yond the scope of this article. This question has however 
been treated in early work on ultrasonic, attenuation in 
metals with anisotropic Fermi surfacesE3. There it was 
found that charge neutrality can be simply imposed in 
terms of an appropriately chosen spatially varying chem- 
ical potential. This leads to the same formula for the 
attenuation as is obtained when the charge-neutrality 
correction is neglected, except that the original electron- 
phonon matrix element is replaced by a related effective 
electron-phonon matrix element. Translated into our no- 
tation, their result is that the electron-phonon matrix 
element Fj(k, q) in Eq. fj] should be replaced by the ef- 
fective electron-phonon matrix element 

F j (k,q)=F j (k,d)-{F j (k,c l )) FS . (9) 

Therefore, in what follows, Eq. with Fj replaced by 
Fj will be assumed to be the correct expression for 
the phonon relaxation rate. It should be noted that, 
for phonons propagating along high-symmetry direc- 
tions, Fj (k, q) differs from Fj (k, q) only for longitudinal 
phonons. 

The phonon lifetime in the normal state can be cal- 
culated from Eq. (Q) by putting the superconducting 
energy gap equal to zero, which gives 



1 



'qj 



8w q,J Tr! 



(if(k,q))__-<F i (k,q)) 



2 

FS 



(10) 



where r„ is the electronJifetime in the normal state. 

In some early workErta on ultrasonic attenuation the 
model Hamiltonian was formulated in terms of the poten- 
tial energy Vipc), assumed to be a function of the contin- 
uous electron position-coordinate x, and describing the 
interaction of the electron with the lattice and with im- 
purities. In this formulation, two problems arise that are 
eliminated by a canonical transformation to a coordinate 
system fixed to the moving lattice. The first is that the 
perturbation of the potential due to the distorted lattice 
is not necessarily small. This problem does not occur in 
the tight-binding formulation because the perturbation is 
naturally formulated in terms of a strain coordinate (i.e. 
the ratio of the displacement to the phonon wavelength) 
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rather than simply a displacement coordinate [see Eq. 
@], The second problem is that it is simpler to work in 
a coordinate system in which the impurities are static in 
the zero-order Hamiltonian so that the electron's energy 
is conserved in the impurity scattering process. In the 
tight-binding approach of this article the impurities are 
static in the zero-order Hamiltonian. For example, poten- 
tial scattering by impurities can be modelled by adding 
an extra term U UiIii cl n . ^c^n^o- to the Hamiltonian for 
the v th orbital on the impurity site labelled by n^. There 
is no dependence of this impurity potential on the ion 
displacements, and the impurities are thus static in this 
representation. Thus, our approach has features quali- 
tatively similar to those of some previous approachesEm, 
even though the details of implementation are different. 

Before proceeding further, we show that our formula- 
tion reduces to the traditional isotropic electron-stress-, 
tensor formulation in an appropriate limit. It is knownL3 
that for hexagonal crystals, sound-wave propagation is 
isotropic with respect to rotations about the c axis. (This 
is not true for tetragonal crystals, such as Sr 2 Ru04.) 
Therefore, we expect that the isotropic electron-stress- 
tensor model could give an appropriate description of 
sound-wave attenuation in a two-dimensional hexagonal 
lattice, particularly if the Fermi surface is taken to be a 
small circle about the point k = so that the electron 
energy spectrum can be approximated by that for free 
electrons, but with an effective mass. For these reasons, 
we consider for the moment a two-dimensional hexagonal 
lattice. The quantity Fj (k, q) defined in Eq. (||) can be 
written in the form 

Fj (k, q) = g X] &* Up (k)e i( g (11) 

a/3 

where a and (3 are summed over the x and y components 
of two-dimensional vectors. Here 

/ Q/3 (k) = J2 RaRfi C0S ( k ' R ) ( 12 ) 
R 

and the sum over R is over three nearest-neighbor hexag- 
onal Bravais-lattice vectors that make angles of 120 de- 
grees with each other. Now, for a Fermi surface that is 
a small circle surrounding the point k = 0, the approxi- 
mation cos(k • R) w 1 — |(k • R) 2 is valid. This leads to 
the result 

/a/3(k) = /a/j(k) - (/ Qj 9(k)>FS 

= -la 2 {k a k p -\k 2 8 a0 ) (13) 

where a is the lattice constant. This is precisely the 
form taken by the effective electron-phonon interaction in 
the isotonic electron-stress-tensor model used by many 
authorsErtil, and shows that our formulation of the the- 
ory of ultrasonic attenuation is equivalent to theirs in an 
appropriate limit. 

To prepare for a more detailed study of the ultrasonic 
attention as observed in Ref. YU explicit expressions for 



Fj(k, q) [derived from Eq. (g)] are now given for our 
Sr 2 Ru04 model for the case that both the phonon wave 
vector q and its polarization vector ej (q) are in the basal 
plane. Furthermore, a simple model of isotropic phonons 
is assumed in which the polarization vector is parallel to 
the wave vector for longitudinal phonons, and perpen- 
dicular to the wave vector for transverse phonons, even 
when the wave vector is not along a high symmetry di- 
rection. The direction of the phonon wave vector q in the 
basal plane is characterized by the angle <j> that it makes 
with the x (i.e. a) axis. 

For longitudinal phonons interacting with electrons in 
the 7 sheet of the Fermi surface via nearest-neighbor in- 
teractions 

F2(k,q) = <7 7 [cos 2 (j)cos(k x a) + sin 2 <pcos(k y a)]. (14) 

For transverse T\ phonons (T% phonons have their po- 
larization as well as their wave vector in the basal plane) 
interacting with electrons in the 7 sheet of the Fermi 
surface via nearest-neighbor interactions 

Fy 1 (k, q) = g 1 cos4>sm(j){cos(k x a) — cos(k y a)]. (15) 

For longitudinal phonons interacting with electrons in 
the xz sheets of the Fermi surface via nearest-neighbor 
interactions 

F[ z (k,q) = g a(3 cos 2 cj)cos(kfa), (16) 

while for interactions with electrons in the yz sheets of 
the Fermi surface 

Ff (k,q) =g^sin 2 <Pco S (kfa). (17) 

In arriving at this result, this hybridization of the xz and 
yz bands has been neglected, as have interband transi- 
tions. The bands are thus one dimensional and charac- 
terized by the Fermi wave vector k^f . 

For transverse T\ phonons interacting with electrons in 
the xz sheets of the Fermi surface via nearest-neighbor 
interactions 

Fff(k,q) = -g aS3 cos sin <j> cos(kf a) , (18) 

while for interactions with electrons in the yz sheets of 
the Fermi surface 

F|*(k,q) = 5 a/3 cos0sin ( /.cos(/c^a). (19) 

Now assume that the attenuation of longitudinal 
phonons is dominated by their interaction with electrons 
on the a and (3 sheets of the Fermi surface as described 
by Eqs. ( [l6|) and ([[7]) (and that the interaction with elec- 
trons of the 7 sheet can be neglected). The attenuation 
will therefore be proportional to 

(Fl) FS - (F L ) 2 FS =p^co S (kf)r 

x (cos 4 + sin 4 cf> - p aP ) (20) 
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Here p a @ is denned to be the fraction of the density of 
states associated with the xz (or the yz) band. (The 
fraction of the density of states associated with the 7 
band is called p 1 , so that p 1 + 2p a/3 = 1). 

It is convenient to define the longitudinal anisotropy to 
be rjioo/rjuo where the mode viscosities of the L100 and 
L110 phonons, tjlwo and 771,110, respectively, are defined 
in the caption to Fig. [j]. The experimental value of the 
longitudinal anisotropy is approximately 30 (see Fig. [I]) . 
The theoretical formula for the longitudinal anisotropy 
in the case that the attenuation is dominated by phonon 
interactions with electrons in the xz and yz bands is [from 
Eq. ©] 
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a/3 



na/3 



(21) 



Using the experimentally determined valuefej p a P = 0.21 
gives a longitudinal anisotropy from Eq. ( [j"l| ) of 2.7. This 
value of the longitudinal anisotropy is too small by a fac- 
tor of 10 to account for the experimentally observed value 
(which is 30). Thus interactions of longitudinal phonons 
with electrons in the a and (3 sheets of the Fermi surface 
can not account for the strong anisotropy of the attenua- 
tion of longitudinal phonons observed in S^RuOJj, and 
these interactions must be dominated by other interac- 
tions. 

Now suppose that the attenuation of longitudinal 
phonons is dominated by their interaction with electrons 
in the 7 band, and neglect the contribution to the at- 
tenuation from electrons in the xz and yz bands. The 
longitudinal anisotropy in this case is given by 



! (/c x a)) 7 - p 1 (cos{k x a))~ 



Vlwo 

Vliio ({i[cos(fc x a) + cos(fc !/ o)]} 2 ) 7 - p"* (cos(k x a))^ 

(22) 

Here the subscript 7 on the angular brackets (as in () 7 ) 
indicates that average is over the 7 sheet of the Fermi sur- 
face only. To begin with, for the purposes of obtaining a 
qualitative understanding of Eq. 22|, neglect (cos(k x a)}-y 
since numerical estimates show that it is small. Now note 
that at the X point of the Fermi surface, k x a = n and 
k y = so the cos(k x a) + cos(k y a) = 0. Because the 
7 sheet of the Fermi surface passes close to this point, 
the value of cos^a) + cos(fc y a) will be small here and 
this will contribute to a large value of the longitudinal 
anisotropy. Furthermore, it should be noted that the X 
point is a saddle point of the electron energy versus k 
surface, so that the Fermi velocity Vu goes to zero there. 
This means that in carrying out the Fermi surface av- 
erage using Eq. (g), points close to the X point will 
be more heavily weighted, thus further enhancing thus 
magnitude of the longitudinal anisotropy. To obtain an 
explicit value for longitudinal anisotropy, we make use 
of the tight-binding approximation to £k given in Ref. 
|l3| , with the parameters given in that article. Thus 
we use, for the electron energy in the 7 band, et = 



E + 2t(cos(k x a) + cos(k y aj) + At' cos(k x a) cos(k y a), with 
parameters (Eg - E F ,t,t') = (-0.4,-0.4,-0.12). This 
allows Eq. ( J2^ ) to be evaluated numerically (including 
now the relatively small effect of a nonzero (cos(k x a)}^), 
giving the result (^lioo/^lho) = 37. This is in rea- 
sonable agreement with the experimentally determined 
value of approximately 30 (see Fig. [j]), considering the 
fact that no attempt was made to adjust the Fermi sur- 
face parameters to improve the agreement, and that the 
inclusion of a relatively small contribution from the less 
anisotropic attenuation due the the xz and yz electrons 
would also reduce the calculated longitudinal anisotropy. 
The conclusion is that electrons on the 7 sheet of the 
Fermi surface give the dominant contribution to the ul- 
trasonic attenuation, and that the detailed Fermi surface 
structure for electrons on the 7 sheet is important for 
understanding the large anisotropy in the attenuation of 
the longitudinal sound waves. 

So far, only the effect of the nearest-neighbor electron- 
phonon interactions on the ultrasonic attenuation has 
been discussed. Note that for these interactions, the 
ultrasonic attenuation of the T100 waves is zero (see 
Eqs. ([[1), (||), and (|9|) for = 0). This is because 
transverse waves propagating in the 100 direction do not 
stretch the nearest-neighbor bonds. Such waves do how- 
ever stretch the next-nearest-neighbor bonds, and atten- 
uation of the T100 waves from the next-nearest-neighbor 
electron-phonon interaction is to be expected. This at- 
tenuation associated with T\ basal-plane phonons inter- 
acting with 7-sheet electrons, can be found in terms of 
the quantity 



^^(k, q) = g'~' cos(20) sin^a) svn(k y a). 



(23) 



The fact that the attenuation of the T100 phonons is 
about 1000 times smaller than that of the L100 or T110 
phonons is a good indication that, in Si'2Ru04 the next- 
nearest-neighbor electron-phonon interactions are in gen- 
eral less important than the near-neighbor interactions. 
The formula 



?7tioo /ff' 7 \ 2 (sin 2 (fc a; o) sin 2 (fc y a)) 



??L100 



FS 



(cos 2 (k x a)) 



(24) 



FS 



(from Eqs. (|14j) and (|23|)), together with the experimen- 
tally measured values of the viscosities ?7tioo an d Vlioo 
can be used to put an upper limit on the magnitude of 
g' 1 1 g 1 of 0.041. This is an upper limit because, as will be 
argued below, a different interaction could be responsible 
for the T100 attenuation. This upper limit is surprisingly 
small, given that the ratio of the next-nearest-neighbor to 
nearest-neighbor hopping matrix elements for the 7 band 
as estimated from the parameters assumed in Ref. |l3| is 
about 0.3. There is, however, no a priori reason why the 
sensitivity of a hopping matrix element to bond stretch- 
ing should be directly proportional to the magnitude of 
the matrix element itself. 

There is at present no experimental information on 
the attenuation of T 2 phonons, i.e. by definition those 
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transverse phonons with their propagation vector q ly- 
ing in the basal plane, and their direction of polariza- 
tion perpendicular to the basal plane. The attenua- 
tion obtained from a consideration of the body-diagonal 
electron-phonon interaction described above is given in 
terms of 

F^ z 2 nz (k, q) = g^ z,rtz [coscf)sm(^k x a) cos(^k y a) 

+ sin (^> cos (| a) sm(hk y a)\ sin(hk z c). (25) 

Some factors of order unity have been absorbed into the 
definition of g£ z ' , ' z . When squared and averaged over the 
Fermi surface, this formula gives an attenuation indepen- 
dent of the direction of q in the basal plane. 

Similarly, the attenuation of T\ phonons by the body- 
diagonal electron-phonon interaction with the xz and yz 
bands is described by 

F^' v *(\l, q) = ag Sz - nz 003(20) cos(±/c z c) 

x sin(|fe x a) sm(^k y a). (26) 

Here, a is a numerical constant of order unity. The in- 
teraction of T\ phonons by the body-diagonal electron- 
phonon interaction with the 7-sheet electrons is described 
by an identical formula, except we put a — 1 and call the 
coupling constant g" 1 . 

IV. ULTRASONIC ATTENUATION IN THE 
SUPERCONDUCTING STATE 

Now that a satisfactory model for the electron-phonon 
interaction has been established, it is possible to pro- 
ceed with confidence to an interpretation of the ultra- 
sonic attenuation measurements^ performed in the su- 
perconducting state of Sr2RuC>4. The basic idea is ta 
use the principles described by Moreno and ColemanEl 
to gain information about the location of the gap nodes 
in Sr2Ru04. Evidence for the existence of gap nodes, 
some of which is presented in Fig. [I] from Ref. [[] has been 
summarized in the introduction. 



a/a 




FIG. 3. Qualitative behavior of the ultrasonic attenuation 
relative to that in the normal state versus temperature rela- 
tive to T c . The attenuation of sound for the case where only 
inactive nodes are present grows more slowly by a factor of 
T 2 relative to that where active nodes are present, at temper- 
atures well below T c . 

At low temperatures in a superconductor with gap 
nodes, only Bogoliubov quasiparticles in the neighbor- 
hood of the nodes are thermally excited, and hence only 
these quasiparticles can interact with the phonons to ab- 
sorb them (and thus to attenuate a sound wave). The 
lifetime of a phonon is given by Eq. (Q). Everything 
depends on the behavior of the electron-phonon matrix 
element Fj (k, q) at the wave vectors k corresponding to 
the nodes. If this matrix element is nonzero at the nodes 
for a particular phonon, then the phonon can interact 
well with the nodal Bogoliubov quasiparticles, and the 
nodes are said to be active for that particular phonon. 
If, on the other hand, the matrix element is zero at the 
nodes for a particular phonon, then the coupling of the 
phonon to Bogoliubov quasiparticle precisely at the node 
is zero, and grows as the distance from the node is in- 
creased. In this case, the nodes are said to be inactive 
for the phonon in question. The ultrasonic attenuation at 
temperatures well below the gap in the case where only 
inactive nodes are present is proportional to T 2 times,the 
ultrasonic attenuation when active nodes are present! 2 ], as 
illustrated in Fig. |[ 

Suppose that the only nodes in the superconducting 
gap are vertical line nodes in (110) planes. Then it is clear 
from Eq. (|l5| ) that the electron-phonon matrix element 
for the T110 phonon is zero at the nodes (since k x = 
±k y for k in a (110) plane), and (110) nodes are thus 
inactive for the T110 phonon. For T100 phonons, the 
nearest-neighbor electron-phonon matrix element of Eq. 
(|15|) is zero for all k because = 0. Thus, for the T100 
phonon, its activity is determined by the next-nearest- 
neighbor interaction given by Eq. (|23|), or Eq. (pq 1 ), 
which shows that the T100 phonon is active for nodes in 
(110) planes. Similarly, Eq. ( |l4| ) can be used to show 
that the electron-phonon matrix element F is nonzero at 
nodes in (110) planes for both LI 00 (0 = 0) and LI 10 
((/> = 7r/4) phonons, showing that both of these modes 
are active. Thus, if the only nodes in the gap are those 
in (110) planes, the attenuation of the T110 sound wave 
would increase significantly more slowly on increasing the 
temperature from zero than that of the T100, L100 or 
LI 10 sound waves. From Fig. |l| this is clearly not the 
case, thus ruling out states that have nodes only in (110) 
planes. In particular this rules out a superconducting 
spin-singlet state of k\ — k y symmetry, as well as a spin- 
triplet /-state with d = z(k x + iky)(k\ — fc 2 ). 

Next consider the case that the only nodes in the su- 
perconducting gap are vertical line nodes in the (100) 
planes. Arguments similar to those given above show 
that these (100) nodes are active for L100, L110, and 
T110 sound waves, but inactive for T100 sound waves. In 
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this case the attenuation of the T100 sound wave should 
increase markedly more slowly with temperature at low 
temperatures than that of the other three sound wave 
types, and this is not the case. This rules out states that 
have only (100) line nodes, such as the /-wave state with 
d — z^kj; ~t~ % ky ^ kg. ky . 

The arguments of the previous paragraphs have ig- 
nored the effects of the a and (3 sheets of the Fermi sur- 
face. Our implementation of our model is inaccurate for 
the behaviour of the electron-phonon matrix element at 
(110) nodes for the a and j3 sheets of the Fermi surface 
because we have not taken into account properly the hy- 
bridization of the xz and yz bands, which is important 
in the [110] directions. A more detailed discussion could 
be given to remedy this deficiency. This will not be nec- 
essary, however, as it will be seen that the more general 
and powerful symmetry arguments of the following sec- 
tion confirm the conclusions reached above. 

It should be noted that the electron-phonon interac- 
tion developed in this article gives results that are signif- 
icantly different from those obtained when the electron- 
phonon interaction is formulated in terms of an isotropic 
electron stress tensor (e.g. see Refs. ||l0 : ll]). This can 
seen, for example, in the top panel of Fig. 2 of Ref. O , 
where a model with gap nodes in (110) planes is ana- 
lyzed, and the L100 mode shows behavior characteristic 
of interaction with inactive nodes, since its calculated at- 
tenuation grows much more slowly as the temperature is 
increased from zero than the calculated attenuation of 
the T100 mode. As noted in Ref. [ll], the intensity of 
the coupling of the L100 phonon to the electrons in the 
isotropic electron-stress-tensor model is proportional to 
the factor (k^, — ^) 2 , which is zero in (110) planes, mak- 
ing the (110) nodes inactive. This inactivity of the (110) 
nodes for the L100 phonon is an "accidental" inactivity, 
i.e. it is a particularity of the isotropic electron-stress- 
tensor electron-phonon interaction which would not be 
present in a more general model and (as demonstrated 
in the follow section) is not required by symmetry. In 
fact, as shown above, the (110) nodes are active for L100 
phonons for the electron-phonon interaction developed 
in this article. Thus it seems clear that the isotropic 
electron-stress-tensor model should not be used in inter- 
pretations of ultrasonic attenuation experiments which 
aim at locating the positions of the gap nodes in an un- 
conventional superconductor. 

Finally note that the electron-phonon matrix element 
for the attenuation of transverse Ti phonons in Sr2Ru04 
(i.e. those that have a propagation vector q in the 
basal plane and a polarization perpendicular to the basal 
plane) is given by Eq. (p5[). The factor sin(fc2c/2) con- 
tained in this expression means that this matrix clement 
is zero for horizonal nodes lying in the plane k z = or 
in on the k z = 2ir/c surface of the Brillouin zone. Thus, 
horizontal nodes at k z = or at k z — 2n/c are inactive 
for Ti phonons, and the comparison of the attenuation 
for T 2 phonons with that for other active phonons can be 



used as a test for such nodes. Unfortunately, it is more 
difficult to find a definitive test for horizontal nodes that 
might appear at k z = tt/c, such as in the proposal of Ref. 
|3(i| , since there does not appear to be any phonon mode 
that is inactive for such nodes. This is related to the fact 
that such nodes are not required by the symmetry of the 
order parameter, but are accidental. 

So far we have investigated only the effects of the rel- 
atively large differences in the temperature dependence 
of the attenuation that are expected to occur as a re- 
sult of the differences between active and inactive nodes. 
Such differences do not appear to occur in the experi- 
mental results currently available for Si^RuOi. There 
is however an intriguing smaller difference in tempera- 
ture dependence noted by Lupien et al.El This is that the 
ratio of the attenuation of the T100 mode in the super- 
conducting state to its normal-state value is somewhat 
larger at low temperatures than that of the other mea- 
sured sound wave modes (see Fig. |l|). This attenuation 
is much smaller than that of the other measured sound 
wave modes, and was accounted for above in terms of a 
second-neighbor intraplanar electron-phonon interaction. 
This is not the only possibility, however. It is possible 
that the interplanar body-diagonal electron-phonon in- 
teraction could be the most important for this mode. 
This interaction has the characteristic feature that the 
squared matrix element contains the factor cos 2 (fc z c/2). 
If there are horizontal nodes close to either the k z — 
or k z = 2t:/c planes where this factor has its maximum, 
then this would give the T100 sound wave modes a boost 
in attenuation at low temperatures relative to the other 
observed nodes, in agreement with experiment. 

In summary, the results of ultrasonic attenuation 
experimental!, when looked at in the light of the above 
theoretical considerations, provide evidence against ver- 
tical line nodes in cither [100] or [110] planes. Also, the 
experiments are consistent with horizontal line nodes. Fi- 
nally, as to the positions of the horizontal line nodes, po- 
sitions close to either the k z = or k z = 2n/c planes 
appear to be favored. 



V. CRYSTALLOGRAPHIC SYMMETRY, THE 
ELECTRON-PHONON INTERACTION, AND 
ACTIVE NODES 

In the last section an explicit model of the electron- 
phonon interaction was used to deduce the activity of 
nodes of different types relative to a given sound wave 
mode. While there is reason to believe that the model 
for the electron-phonon interaction is a relatively good 
one, since it accounts well for a number of unusual fea- 
tures of the sound attenuation in the normal state of 
Sr2RuC>4, it is nevertheless of interest to develop an ap- 
proach that is not dependent on the details of the partic- 
ular model. Such an approach, which exploits the crys- 
tallographic symmetry, is indeed possible, and will now 
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be sketched. Such symmetry arguments will be particu- 
larly useful in cases where the electronic structure may 
not be well known, or is sufficiently complicated that a 
detailed model is difficult to develop. 

In this section, the Hamiltonian describing the 
electron-phonon interaction will be taken to be that of 
Eq. (^), where the matrix element 5^1 ^ will be as- 
sumed to be completely general, and subject only to the 
restrictions imposed by symmetry. For example, from 
the properties of the Hermiticity of the Hamiltonian, its 
time-reversal invariance, and its invariance with respect 
to spatial inversion in a ruthenium ion position, it follows 
that <?£ + i q q j is pure imaginary, is an even function of its 

argument k+ ^q, and is an odd function of its argument 

q 

Symmetry arguments can be used to show that, for 
a transverse phonon with wave vector q along the [110] 
direction, and its polarization j in the basal plane, per- 
pendicular to q, vertical lines nodes in (110) planes are 
inactive. Suppose that k has its basal-plane component 
in the [110] direction. Then under a reflection in a plane 
normal to the basal plane, and containing q, c v ^ must 
have a definite parity (i.e. either c v ^ a — ► +c u ^ a or 
c ^,k,o- — > ~Cv^a). Also, Aqj — * — Aj. j ■ From these 
properties it follows that 5^ + i q . = for all k lying in 

the same (110) plane that contains q. A related argu- 
ment shows that for k lying in the plane perpendicular 
to q, 5^ qj = 0. For this latter argument, it is neces- 
sary to assume that |q| <C |k|, so that can ^ e 
replaced by g£qj- Thus, for a T110 sound wave, verti- 
cal line nodes in (110) planes are inactive. Furthermore, 
there are no symmetry arguments that show that (110) 
nodes are inactive for T100, L100 or L110 sound wave, 
and they must therefore, in general, be active. 

In a similar way it can be shown that symmetry re- 
quires that (100) vertical nodes are inactive for T100 
sound waves, while these nodes are in general active for 
T110, L100 and L110 sound waves. 

Making use of the fact that the basal plane is a plane 
of reflection symmetry, one can also show that line nodes 
in the plane k z = 0, and in the plane k z — 2?t/c are 
inactive for Ti sound waves having q in the basal plane 
and polarized perpendicular to the basal plane. 

In the long- wavelength limit |q| <C |k|, the symmetry 
arguments given above can be formulated in a particu- 
larly simple form. In this limit, the function _FJ(k, q) 
determining the symmetry of the electron-phonon inter- 
action [and which can be obtained from from Eqs. (^) 
and (^)]can be written as 

F/(k,q)= ff " J2 fZfiMteM- (27) 

cn,/3— x,y ,z 

(Note that, in the isotropic electron-stress-tensor model, 
faff fe) — k a kp — (l/d)6 a /3, where d is the dimensionality 
of the system; this isotropic stress-tensor expression is 
not used here.) If, for a given k, there is a symmetry 



operation from the crystallographic point group which 
leaves k invariant, but changes sign of q a e.j,@, then the 
matrix element /^(k) vanishes (one should also use the 
fact that in symmetry operations, Bj and q transform 
like vectors, and Bj(— q) = ej(q)). In particular, one of 
the consequences is that there are no zeros of i<J(k, q) 
required by symmetry for all longitudinal waves, so the 
nodes are always active for this polarization. Also, it is 
clear that a symmetry-imposed line zeros for transverse 
waves can be present only if k is in a high-symmetry 
plane. 

Interband transitions of electronic quasiparticles have 
been ignored in our discussion. Such transitions are ex- 
pected to play a role only when two bands cross, and 
since there are relatively few quasiparticles associated 
with such points, their effects in general should not be 
important. However, if a case is encountered where inter- 
band transitions play a significant role they will have to 
be investigated carefully, since the symmetry properties 
of the electron-phonon matrix element for such transi- 
tions are different from those for intraband transitions. 

As an application of these ideas to. *u material with a 
relatively complicated Fermi^rfaceEj'LJ, consider ultra- 
sonic attenuation in UPt3,OCj which is believed to give 
evidence for the existence of basal-plane line nodes in the 
superconducting gap. We examine the consequences of 
the conjecture that the order parameter is characterized 
by a singlet-state gap A(k x , k y , k z ) transforming in ac- 
cordance with the Fi2 g representation of the point group 
Dg/j, or by the triplet-state order parameter d z {k Xl k y ,k z ) 
which transforms according to E2 U . (These are two of 
the most commonly assumed candidates for the order 
parameter, e.g. see Ref. [n].) First note that from their 
symmetry classification, both of these order parameters 
change sign under reflection in the basal plane. (The 
space group of UPt3 is P63/mmc, which has a a z reflec- 
tion plane.) Hence, 

A(fc ;E , ky, k z ) — A(k Xl ky, k z ) 

d z (k x , ky , k z ) = d z (k x , ky, fc z ) . (28) 

The constraints of the this equation require the order 
parameters to be zero at k z — 0, i.e. there will be a 
line of nodes wherever a sheet of the Fermi surface cuts 
the basaLplane. (The Fermi surface consists of several 
sheetsES'Ej, with each sheet having its own gap, but all 
gaps are expected to have the same symmetry, and to sat- 
isfy Eq. (p28|).) Note that Eq. ( ^8|) also requires the order 
parameters to be zero at surfaces of the hexagonal Bril- 
louin zone at k z = ±7r/c. This is because {k x , k y , +7r / c) 
and (k x , k y , ~tt / c) are equivalent points. There are two 
toroidal sheets of the Fermi surface, that cut the surface 
k z = tt/c of the Brillouin zone P 8 " 39 ! and hence give line 
nodes there. The presence of these Brillouin-zone sur- 
face line nodes, associated by symmetry with those in 
the basal plane, does not seem to have been previously 
noticed. 

Because the basal plane of UPt3 is a plane of reflec- 
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tion symmetry, symmetry arguments similar to those 
given above show that both basal-plane line nodes and 
k z = tt/c Brillouin-zone-surface line nodes are inactive 
for T2 sound waves with q || a and polarization parallel 
to c, but active for T\ sound waves with q || a and po- 
larization parallel to b. The tem pera ture dependence 
of the ultrasonic attenuation datac2lei shows that the 
nodes are active for T\ sound waves, and inactive for 
T 2 sound waves, giving clear evidence of the existence 
of basal plane and k z = n/c Brillouin-zone-surface line 
nodes. These comments extend previous discussions to 
include Brillouin-zone-surface line nodes, and also ensure 
that the analysis remains valid for the complex Fermi sur- 
face geometry that actually occurs in UPt3. 

VI. CONCLUSIONS 



method is useful even when it is not possible to develop 
a detailed microscopic model for the electron-phonon in- 
teraction. 

As a by-product of our theory, we also propose that the 
ultrasonic attenuation data in UPt3 can be interpreted 
in favor of the existence of horizontal line nodes in the 
plane k z = ±7r/c, as well as in the plane k z = 0. 

Our results, both for the anisotropy of the ultrasonic 
attenuation in the normal state, and for questions re- 
lated to the positions of gap nodes in the superconduct- 
ing state, are different from those obtained when using 
the traditional method of modeling the electron-phonon 
interaction is terms of an isotropic electron stress tensor. 
In particular, the isotropic electron-stress-tensor model 
makes certain superconducting gap nodes accidentally in- 
active, when they are in fact active. 



The electron-phonon interaction in Sr 2 Ru04 has a 
strong anisotropy that is highly unusual. The unusual 
nature of this interaction can be attributed to certain 
structural properties of Sr2Ru04, these being that the 
structure is a layered one with a relatively large dis- 
tance between the RUO4 layers, and that the Ru ions 
in a layer form a square lattice with the largest electron- 
phonon interaction being between nearest-neighbor Ru 
ions. A detailed model electron-phonon interaction based 
on the idea of a tight-binding Hamiltonian with hop- 
ping matrix elements that depend on the ion displace- 
ments does a good job of quantitatively accounting for 
the huge anisotropy observed^ in ultrasonic attenuation 
in the normal state of Sr2Ru04. The dominant con- 
tribution to the attenuation comes from the interac- 
tion of phonons with electrons in the 7 band. The 
attenuation of transverse sound waves propagating in 
the [100] direction and having their polarization in the 
basal plane is exceptionally small because such waves 
do not stretch the nearest-neighbor bond lengths, and 
hence have no nearest-neighbor electron-phonon interac- 
tion. The strong anisotropy for the propagation of lon- 
gitudinal waves in the basal plane can be related to the 
fact that the 7 sheet of the two-dimensional Fermi sur- 
face describes a large circle passing close to the X point 
of the Brillouia zone. An analysis of the ultrasonic at- 
tenuation datatl in the superconducting state in terms of 
our model electron-phonon interaction rules out the pos- 
sibilities that the only nodes in the superconducting gap 
are vertical lines nodes in (100) planes, or vertical line 
nodes in (110) planes. The experiments performed so far 
are consistent with the existence of horizontal line nodes 
such as those proposed in Refs. and 30. With respect 



to the positions of horizontal line nodes, positions close 
to the planes k z = or k z = 2tt/c appear to be favored, 
but a more detailed analysis is need to confirm this. 

A general method (based on crystallographic symme- 
try arguments) of determining the existence of inactive 
nodes for a given sound wave mode is developed. This 
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APPENDIX: SYMMETRY AND UNIVERSALITY 
OF ULTRASONIC ATTENUATION 

One of the remarkable features of unconventional su- 
perconductors is that some of their transport coefficients 
attain "universal" (i.e. independent of the disorder con- 
centration) values at sufficiently low temperatures. In 
particular, universality was first predicted for electrical 
conductivity of d-wave superconductors in Ref. WA ( sec 
also Ref. [43|) , and for thermal conductivity in Refs. |44|j45| 
and observed experimentally in thermal conductivity in 
Ref. In this Appendix, we study the low-temperature 
behaviour of the ultrasonic attenuation in unconventional 
superconductors, using the model formulated in Sec. [n]. 

In the absence of the vertex corrections, the ultrasonic 
attenuation coefficient in the hydrodynamic approxima- 
tion is given by 

aj (q,T c ) r n J Q a£ { de) Ree(e)' 
where j = L, T\, T 2 labels the phonon polarization, and 
1 



x(F?(k,q)Re 



2Ime(e) \ J ' ' " / fs 

' g( e )| 2 + g 2 ( £ )-2|A k | 2 \ 



(30) 
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The temperature-independent ultrasonic attenuation in 
the normal state is given by 



aj(q,T c ) 



8w 



N F T n F/(k,q) 



FS 



(31) 



The expression (^9|) is similar to that given in Ref. | the 
only difference being in the angular dependence of the 
electron-phonon interaction (see also Ref. |]). To repro- 
duce the results obtained with the help of the isotropic 
electron stress tensor, one should replace Fj (k, q) — > 
(k.q)(k.e 3 (q)) 7 (l/d)(q-e,(q)). 

The function e(e) describes the self-consistent renor- 
malization of quasiparticle energy due to impurity scat- 
tering and satisfies the equation 



_J #(£) 

2t„ cos 2 S + g 2 (e) sin 2 6 ' 



where g(e) is the retarded Green's function G at coin- 
ciding points (see Ref. |l8| for a review of the effects of 
disorder in unconventional superconductors). Assuming 
electron-hole symmetry, 



9(e) 



v^ 2 -|A k | 2 



(33) 



In particular, for quasi-2D d-wave order parameters 
Ak = Aocos2<ys (d x 2_ y 2 symmetry) or Ak = Aosin2<^ 
{d X y symmetry) and a cylindrical Fermi surface, g(e) = 
(2/n)K(Ao/e), where K(x) is the complete elliptic in- 
tegral. The function e(e) determines, for example, the 
inverse mean free time of quasiparticles: 



2Ime(e), 



(34) 



and the disorder-averaged quasiparticle density of states 
(DoS): 



N(e) = N F Reg(e(e)). 



(35) 



The behaviour of the attenuation coefficient strongly 
depends on the value of the phase shift So . Here we con- 
sider two limiting cases of weak impurities (Born limit), 
when Sq — > 0, and strong impurities (unitary limit), when 
S a -> vr/2. 

The solution of Eq. ( |32|) at zero energy is purely imag- 
inary in both cases: e(e = 0) = iTo, where 



T = A cxp(-7rA r„) 
in the Born limit, and 



r = A 



2A t„ In A t„ 



(36) 



(37) 



in the unitary limit. The zero-energy scattering rate 
Tq determines the crossover energy scale separating two 



qualitatively different types of the behaviour of the ob- 
servable quantities. If the typical energy of excitations 
(temperature) is greater than To, then one can neglect 
the self-consistent energy renormalization and use the 
quasiparticle Boltzmann equation for calculating the ki- 
netic properties in the superconducting state (for the ap- 
plication of this approach to unconventional supercon- 
ductors, see, e.g. Ref. In contrast, if the typical 
energy is smaller than To, then the self-consistency ef- 
fects become important, and the quasiparticle Boltzmann 
equation is not applicable. In the former case, the imag- 
inary part of e is small compared to Re e — ► e, and we 
obtain 



,-(q,r) 



i 



Oj(q,T c ) r„ 



de - 



df 
de 



M<1> e ) 



(32) where 



A J -(q,e) = 2r(e) 



(i? 2 (k,q)R Cv /e 2 -|A k || 2 ) 

\ / F S 



(38) 



(39) 



The expression © is equivalent to Eq. (0) for isotropic 
impurity scattering. It should be noted that the quasi- 
particle scattering rate ( |34|) in the Born limit decreases 
with energy in the superconducting state, and at some 
temperature T* ~ (wr„)T c the applicability condition of 
the hydrodynamic approximation is violated. It turns 
out, however, that in real experimental conditions for 
Sr2Ru04, the crossover temperature T* is so small that 
we neglect this complication here. 

In the low-temperature regime, we replace Ime by Tq, 
take the limit Ree — + 0, calculate the integral over en- 
ergy, and obtain 



<*j(q,r) = J_ 
a,(q,T c ) 2r n ' 



^(k,q) 



r 2 
1 



(r 2 + |A k | 2 ) 3 / 2 



FS 



\ J I FS 



(40) 



We see that, at T < To, the ultrasonic attenuation does 
not depend on temperature. 

In the Born limit, for typical disorder concentrations 
in cuprates and ruthenates, F turns out to be exponen- 
tially small compared to T c (in particular, in the exper- 
imental conditions of Ref. |l|, Aot„ ~ l n .ab/£,o,ab ~ 10, 
so that Tq/T c < 1CP 3 ). For this reason, only the "high- 
temperature" limit is relevant for Born impurities. 

For unitary impurities, however, Tq can be as large as 
0.1T C , which means that, on one hand, there might be no 
clear power- law behaviour of a(T) at Tq < T <C T c (it 
was pointed out, e.g. in Ref. ||), and, on the other hand, 
the low-temperature regime can be observable. The con- 
tributions to the ultrasonic attenuation from active and 
inactive nodes can be easily separated because of their 
different dependences on the impurity concentration. For 
active line nodes at a cylindrical Fermi surface, 
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if (k,q) 



r 2 
1 o 



F 2 



(r 2 + |A k p)3/ 2j 



FS 



1 

7T Ao ' 



(41) 



where Fo is the value of Fj (k, q) at the line of nodes. 
Therefore, as seen from Eqs. ( |3l| ) and (40), the ultra- 
sonic attenuation does not depend on the impurity con- 
centration. In contrast, for inactive line nodes, 



if(k,q) 



r 2 

1 



(r 2 + |A k | 2 )3/2 



FS 



A 3 - 



7T 



*0 







(42) 



where Fq is the value of the transverse derivative of 
Fj(k, q) at the line of nodes. Therefore, in this case, 
the ultrasonic attenuation does de pen d on the impurity 
concentration. Comparing Eqs. ( ]IT|) and (ft^), we see 
that the contribution from inactive nodes is typically 
much smaller than that from active ones, and one can 
make a conclusion that if, for a given polarization and 
propagation, there are active nodes present, the atten- 
uation coefficient is universal at low temperatures. In 
particular, the attenuation of longitudinal waves should 
always be universal. On the other hand, for example, 
the attenuation of the in-plane T2 waves for the order 
parameters with horizontal line nodes at k z — cannot 
be universal, because such nodes are inactive. Also, for 
a /-wave order parameter d(k) oc i{k x + ik y )k x k y , the 
attenuation coefficient of the T110 phonons is universal, 
whereas that of the T100 phonons is not. For the order 
parameter d(k) oc z(k x + ik y )(k x — fc 2 ) the situation is 
opposite: the T100 attenuation is universal, but T110 is 
not. It should be noted that the calculations based on 
the isotropic electron-stress-tensor model give different 
results! 
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